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1. Introduction 

Sphaleron processes are important for baryogenesis PJ El El Ej This is obviously true in 
the case of electroweak baryogenesis because the weak sphalerons are the only processes 
that violate baryon number. But is also the case of leptogenesis mechanisms 00013, 
where weak sphalerons are responsible for converting part of the generated B — L 
asymmetry into a baryon asymmetry. In this case the effect of weak sphalerons can be 
reduced to a mere constraint on a linear combination t] w of quark and lepton chemical 
potentials [HI GUI E] of the form r] w = 0. But the constraint only applies when sphalerons 
are in thermal equilibrium in a range of temperatures from around 10 12 GeV to the 
electroweak phase transition. It is important to know well this temperature range when 
designing a lepto-baryogenesis mechanism, as well as all the other reactions in thermal 
equilibrium at every moment as the Universe cools down, because too many constraints 
enforce a trivial particle asymmetry result. On the other hand, the surviving baryon 
asymmetry depends quantitatively on the consecutive sets of constraints applying at 
and after the epoch of baryogenesis or leptogenesis. 

The weak sphaleron frequency equalizes the Hubble expansion rate at a certain 
temperature. However, it cannot be identified with the relaxation temperature of the 
sphaleron chemical potential constraint. It is the main purpose of this paper to study 
the evolution of r/ w and evaluate its relaxation temperature as well as for the QCD 
sphaleron ^21 EH! degeneracy parameter, r] B . It is also argued that the running of the 
weak and strong coupling constants with the energy scale leads to smaller values of 
the sphaleron rates at very high energies, specially in the case of strong sphalerons. 
In the next section we establish the general relation between charge violation rates 
and sphaleron diffusion rates using standard scattering theory considerations. In the 
subsequent sections we specialize to the cases of weak and strong sphalerons obtaining 
general results for the damping rates of the parameters r/ w and r] s . We calculate also 
the respective relaxation temperatures in the standard model case and in a typical 
leptogenesis scenario. In section H] we analyze the interplay between weak sphalerons, 
bottom quark an tau lepton Yukawa couplings while they approach thermal equilibrium 
and derive the appropriate system of coupled equations. In section |H1 we extend the 
study to the context of supersymmetric standard model versions and in the last section 
summarize the conclusions. 

2. Sphaleron processes 

The rate per unit of volume of any reaction with one or more particles a, in the initial 
state and bj in the final state, a, — > bj, depends on the transition amplitude T and 
particle distribution functions f a as follows [T^ ITo] Ho]: 




(1) 
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where a product over i and j is implicit. <i$ stands for the phase space element and 
1 ± are the final state boson stimulated emission and/or fermion Pauli blocking 
factors. The evolution of any charge asymmetry depends on the difference between 
the rates of a, — > bj and its inverse reaction bj — > dj, 7+ and 7_ respectively. Using 
Bose-Einstein and Fermi-Dirac distribution functions and assuming that the transition 
amplitude is invariant under time reversal, one obtains to first order in the chemical 
potentials, 

7 + - 7- = ~(7 + +7-)£ Y ANa ' (2) 

a 

where AN a is +1 (—1) for each of the outgoing (incoming) particles. An observable Q 
that varies AQ = Q a AN a (-AQ) in a single reaction a, — > bj (bj —>■ Oj) suffers a net 
violation per unit of time in a total volume V equal to dQ/dt = (7+ — 7_) VAQ. 

Sphalerons processes produce or destroy fermion particles and are characterized 
by a Chern-Simons number increase AiV whose sign changes with the direction of the 
reaction as any other charge. Then, the violation rate of any charge Q is obtained by 
summing over all AiV numbers as follows: 



dQ 



|£ 7 AivAQ£^AiV a T/. (3) 



dt 2 f-f '" JV ^— ' T 

AN a 

It happens that the fermion number variations AN a are directly proportional to the 
Chern-Simons number variation AiV, so that AQ/ AN and the degeneracy parameter 

v = y^L^ (4) 

^ T AN K ' 

a 

are independent of AiV. This gives 
dQ 1AQ 

-dt=-2AN TvV > (5) 
r = £ 7 A7v(AiV) 2 , (6) 

AN 

where T is nothing but the sphaleron Chern-Simons number diffusion rate calculated in 
the literature 0E7], 

r= Urn mt)-N(0m 

V,t->oo Vt V ' 

3. Weak sphalerons 

Weak sphaleron processes of Chern-Simons number variation AiV produce fermion 
number violations [T5J HHI equal to AiV^ = AiV per lepton doublet k and AN qi = 
AN per quark doublet and colour state where % denotes the fermion generation. 
Specializing the general expressions (@J) and (0) to weak sphalerons, one derives the 
baryon and lepton number violation rates per generation, 
dB a dLi. 1 
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where the degeneracy parameter, 

r? w = ^^(3// 3i + /^.) , ( 9 ) 

i 

is a function of the quark and lepton doublet chemical potentials, degenerate in colour 
and isospin during the electroweak symmetric phase. One recovers so the known result 
for the baryon number violation rate [01 E0J EQ (notice that |r w represents the average 
between the forward (AiV > 0) and backward sphaleron rates). 

The detailed balance equation r] w = holds if weak sphalerons are in thermal 
equilibrium, but not otherwise. rj w is a dynamical damped quantity whose relaxation 
time is much larger than the age of the Universe when sphalerons are not in equilibrium 
and much smaller when they are. The r] w evolution depends not only on sphalerons 
but also on all the other reactions occurring at that moment. The rapid reactions 
enforce detailed balance equations that reduce the number of independent chemical 
potentials. The latter have a one-to-one relation with an equal number of independent 
quantum numbers that are either conserved or are violated solely by reactions not fast 
enough to be in full thermal equilibrium. Without loss of generality, one can choose 
the independent quantum numbers as the baryon number B and a set of charges, Q a , 
that are conserved by weak sphalerons: if a given quantum number is violated by weak 
sphalerons there is always a linear combination with B, for example B — L, that is 
conserved by them. r] w can be written in terms of B and Q a as 

r] w T 3 V = aB + b a Q a . (10) 

After reheating the total entropy is conserved and so is the product T 3 V, whenever the 
number of relativistic degrees of freedom is stable. Under these conditions, and taking 
the usual approximate expressions, linear in the chemical potentials, of the particle anti- 
particle number asymmetries, the coefficients a and b a are well defined constants that 
only depend on the kind of particles and reactions in thermal equilibrium. They are 
thus model dependent but do not depend on the origin of the particle asymmetries. 

Let us assume that at the temperature range of interest weak sphalerons are 
the only processes that violate baryon number and that all reactions besides them 
can be classified as either rapid or negligible slow. The first enforce the chemical 
potential constraints, the latter do not contribute to flavour or any other quantum 
number violation rates. In this case the charges Q a are exactly conserved quantities 
by construction, whose values are determined by the initial conditions and particular 
baryogenesis or leptogenesis mechanism. Then, one derives from equations (jHJ) and (jTUj) 
that for 3 fermion generations, 

d,Tj w 3 T w . . 

The sphaleron rate goes as T w ~ T 4 and in a radiation dominated Universe the Hubble 
expansion rate behaves as H = l/2t ~ T 2 . Hence, T w = T w T~ 2 H~ l is a constant 
quantity. T w may be naively considered as the sphaleron equilibrium temperature 
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because T W T 3 = H at T = T w , but this is an oversimplification. In fact, one obtains 
that i] w has a relaxation temperature equal to T* = 3/2 a T w : 

^pi = -^aT w T] w . (12) 

The parameter t/ w decays exponentially as e~ T ™/ T and its initial value follows from 
equation (jl(J|) with the initial baryon number asymmetry. As r/ w decays, B converges to 
the value that annihilates the right-hand member of equation (|lUj) . 

The coefficients a and b a are model dependent. They depend on the particle content 
and reactions in equilibrium. The standard model interactions that enter in equilibrium 
before the 10 12 GeV temperature scale are the gauge interactions, implying zero gauge 
boson chemical potentials and fermion chemical potentials degenerate in colour and 
isospin, QCD sphalerons O EH], an d the top quark Yukawa coupling (q t t c R (j)). They 
put the following constraints on the Higgs doublet 0, up and down quark isosinglets 
Ui = ujt,CR,tn, di = dR,SR,bR, and quark doublets qi = q u ,q c ,Qt defined as having 
diagonal Yukawa couplings with the up quarks: 

5^(2/i 9l - ~ MdJ = , (13) 

i 

/% - Vt R + H = , (14) 
The other constraint is the initial condition of zero hypercharge, 

^2(p qi + 2/i Ui - iL d . - fi h - fi e J + 2^ = , (15) 

i 

to first order in the chemical potentials and neglecting thermal masses. Above the 
10 12 GeV temperature scale the other quark and lepton Yukawa couplings are negligible 
and the gauge interactions, including sphaleron processes, conserve several flavour 
quantum numbers, namely: lepton numbers L ei of e« = 6r, /jlr,tr; differences between 
the baryon numbers of right-handed charm or down quarks and the Ur baryon number, 
B CR — B UR and B^ — B UR ; differences between the total lepton or baryon numbers 
of different generations, — Lj and Bj — Bj; B — L. While B — L and Lj — Lj are 
conserved by all standard model interactions the other quantum numbers are violated by 
Yukawa couplings. The up quark is singled out because it has the smallest quark Yukawa 
coupling. As the Universe cools down and the smaller Yukawas enter in equilibrium the 
various flavour charges cease to be conserved and get constrained by the new detailed 
balance equations. 

If there are non-standard interactions in equilibrium, such as lepton number 
violating heavy neutrino Majorana masses, more constraints apply and less independent 
charges subsist. The remaining conserved charges at any given moment, Q a , are 
determined by the particular baryogenesis or leptogenesis mechanism, but their values 
do not affect the r] w damping rate (not even the zero hypercharge condition). It only 
depends on the constraints enforced by the interactions in full equilibrium. Assuming 
that only the standard model gauge interactions, strong sphalerons and top quark 
Yukawa coupling are in thermal equilibrium, one obtains a = 15/2 for the coefficient a 
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in equation JTUj). It follows that the r] w relaxation temperature is T* = 45/4 T w , one 
order of magnitude larger than T w . 

The sphaleron diffusion rate has been shown to depend on the temperature 
and weak coupling constant as T w « 26 o^T A . At this point we argue that the effective 
coupling constant runs with the energy scale and therefore varies with the temperature \. 
At the 10 12 GeV energy scale a w « 1/41 which gives r w « 2 x 10~ 7 T 4 , about 5 
times smaller than the rate calculated with the Mz scale value a w ~ 1/30. For a 
Hubble rate H = 1.66 g l J 2 T 2 / M?\ with m 107 effective number of degrees of freedom 
the T w and relaxation temperature values are respectively T w « 1.6 x 10 11 GeV and 
T* » 1.8 x 10 12 GeV. 

Leptogenesis is a mechanism that has been actively investigated IE] • It is 
based on the addition of extra sterile neutrinos iV a with very heavy Majorana masses and 
complex Yukawa couplings, kN a (j), with the standard Higgs and lepton doublets. They 
cause lepton flavour and total lepton number violating processes that stay in thermal 
equilibrium for a certain period of time. During that period the following constraints 
apply: 

#4+^ = 0. (16) 

The quantum numbers Lj — Lj and B — L are then rapidly violated and consequently 
they are excluded from the set of charges Q a in equation (JTUJ). After applying these 
new constraints one obtains for the coefficient a the value a = 261/44 and the slightly 
reduced relaxation temperature T* ~ 1.4 x 10 12 GeV. 



4. bji and tr Yukawa couplings and sphalerons 

The weak sphaleron relaxation temperature falls in the vicinity of the equilibrium 
temperatures of the and tr Yukawa couplings. Our estimates of the bn and tr 
Yukawa production rates Y b and T T i.e., reactions where &r and r R are produced in 
the final state, yield equilibrium temperatures equal to = T&T H^ 1 2 x 10 12 GeV 
and T T = T T TH~ l m 10 12 GeV. This shows that one cannot neglect the violation of 
Q b = Bh R — B UR and L TR when weak sphalerons approach thermal equilibrium. The 
violation rates obtained from the basic expression (j2J are given by 

^± = -2T b n eqVYb V , (17) 
dL 



-^L = -2T T n eqVYT V , (18) 

where n eq = 0.90T 3 /7r 2 denotes the ultra relativistic fermion equilibrium density and 
the degeneracy parameters are defined as 

VY b = Vb R + ~ Vqt , ( 19 ) 

Vy t = Vr R + V<t> ~ Vir ■ (20) 
I I thank M. Shaposhnikov and D. Bodeker for discussions on this point. 
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This goes in parallel with the baryon number violation rate, 
dB 



It 

dt 2 v 



VwV 



(21) 



assuming that B is only significantly violated by sphalerons. 
The above equations have the general form, 
dQf, 



dt 



Vim v 



(22) 



where Q M = (B,Q b ,L TR ) and rj^ = (Vw,VY b ,VY T )- The charges and the set 
of conserved quantum numbers Q a have a one-to-one relation with the independent 
chemical potentials left unconstrained by the detailed balance equations. One can thus 
write the degeneracy parameters 77^ as linear combinations of and Q a : 

% T 3 V = a»Q„ + b«Q a . (23) 

The differentiation against time yields the general result 



-dT = - a ^ T 



Vv 



In a radiation dominated Universe this coupled system of equations leads to: 



dT~ 

rpv 



(24) 

(25) 
(26) 



where is the relaxation temperature matrix. 

If all interactions in equilibrium at around 10 12 GeV are standard model 



interactions, the conserved charges Q a entering in equation (|z3jl are: L 



B Cr B Ur , 



B SR — B UR , Bd R — B UR and Bi — Bj, only negligibly violated by Yukawa 
couplings, and the exactly conserved numbers Lj — Lj and B — L. The constraints 
imposed by the tR Yukawa processes, strong sphalerons and gauge interactions determine 
the coefficients and b". In particular, 



V 



^ -i \ 

w 23 23 

1 12 38 

" 6 23 23 



(27) 



The estimates of the reaction rates in equations ([T7jl - (|22)l given above and at the end 
of the previous section produce the following relaxation temperature matrix: 



rj-ip 



( 



\ 







4 i 

20 



2 




x 10 12 GeV . 



(2f 



The elements 11 of these two matrices coincide with the values of the coefficient a 
in equation (fTUj) and the sphaleron relaxation temperature calculated in the previous 
section, where the 6r and tr Yukawa couplings were neglected. Now, r?^ = r/ w , 
rjy b , rjY T obey the coupled system of equations (J2"5J) . It can be integrated in closed 
form: diagonalizing with a non-orthogonal matrix V, the general solution is, in 
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matrix notation, rj = Ve T *( T ® ~ T 1 ^V~ 1 r]o, where T* is the diagonalized relaxation 
temperature. The rj^ parameters oscillate into each other as they decay but the 
outcome of this complex evolution is clear: weak sphalerons, or and tr Yukawa 
couplings enter in equilibrium at the same time with relaxation temperatures in the 
range (2 - 4) x 10 12 GeV. 

This exercise can be repeated for periods characteristic of leptogenesis in which 
lepton number is rapidly violated by heavy neutrino Majorana masses. Then, the 
additional constraints (fTT)|) are enforced and the conserved charges Q a in equation ()23j) 

By Now, the results 



reduce to: L ea , L m , B CR 



B UR , B Sr 



Bur i B t i R 



B UR and B { 



for the coefficients a v and relaxation temperature matrix are 



3 

22 



/ 



87 
2 

-4 



-36 
64 
24 



-21 \ 



(29) 



rpv 
A* 



V 



-2 
3 
1 



i 

5 
3 
2 



x 10 12 GeV 



(30) 



in the same range of temperatures as the pure standard model case. 



5. QCD sphalerons 

The value of the QCD sphaleron diffusion rate, T s , is still uncertain but there are 
indications (221123! that it scales with the weak sphaleron rate as T s ~ 12«g/a^r w . At 
high energies the strong and weak coupling constants converge to each other: a s ~ 1/37 
and a w « 1/42 at 10 13 GeV, which gives T s ~ 4 x 10~ 6 T 4 , only 20 times larger than T w . 
The reference temperature in a radiation dominated Universe, T s = r s T~ 2 if -1 , takes 
the value T s ~ 3 x 10 12 GeV. Again, this is not the relaxation temperature, T*, of the 
strong sphaleron degeneracy parameter, 

Vs = Tj, ^2( 2 ^ - Vui ~ Vdi) ■ (31) 

i 

In order to study the r] s dynamical evolution one needs the rates of chiral charge 
violation induced by QCD sphalerons. A sphaleron with a Chern-Simons number 
increase AN induces a chiral charge variation AQs = 2AiV for each of the 6 Dirac 
quarks fTHl I19 | IT2" | I13j. defining chiral charge as +1 for left-handed and —1 for right- 
handed quarks, which means a quark number increase AN a = ±AiV/3 per colour and 
chiral state. Employing the general result of equations (jlj), (JSJ) and the violation 
rates of left-handed (qi) and right-handed (qn) quark number asymmetries are 

dt dt 2 ' 1 v ' 

per isospin state (summing over colours). This is in agreement with ref. [22J (again, |r a 
represents the average between forward and backward sphaleron rates). 
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Neglecting thermal mass corrections, rj s is directly related with the total chiral 
charge as Q5 = ^r] s T 3 V, but that is not enough to determine the rj s evolution from 
equation (p?2*|) because the third generation chiral charge is rapidly violated by the top 
quark Yukawa coupling. It is more convenient to express i] s in terms of the right-handed 
up quark baryon number B UR = N UR /3, for example (ur is the quark with smallest 
Yukawa coupling), and a set of independent quantum numbers Q a conserved by strong 
sphalerons: the baryon number B and other charges conserved by both strong and weak 
sphalerons: 

Vs T 3 V = -aB UR + b a Q a . (33) 

The coefficients a, b a and the set of charges Q a depend on which interactions are in 
thermal equilibrium. However, if all processes that are not in full thermal equilibrium 
including weak sphalerons can be considered negligibly slow, the charges Q a are 
conserved in first approximation and only B UR is violated in the right-hand side of 
equation (J3llj) . One derives 

d Vs a r s 

¥ = l^ s ' (34) 

which implies a relaxation temperature equal to T* = aT s /6 in a radiation dominated 
Universe: 

£r = ->. (35) 

If no interactions beyond standard model are in thermal equilibrium, the chemical 
potential constraints are determined by the gauge interactions, top quark Yukawa 
coupling and the condition Y = (equations (JUjl and (jl5|)) . The approximately or 
exactly conserved charges Q a in equation (J31!)) are B, B — L, L e ., B CR — B UR , B d . — B UR , 
Li — Lj and Bi — Bj. In that case one calculates a = 184/3 and the consequent relaxation 
temperature is T* = 10.2 T s ~ 3 x 10 13 GeV. It compares with the value obtained if we 
had ignored the top quark Yukawa coupling, T* = 12 T s . 

% decays as e~ T ^ T and its initial value is obtained from equation (|3^jl with the 
initial B UR value. As 7] s decays, B UR converges to the value that annihilates the 
right-hand side of equation ()33j) . When the weak sphalerons enter in equilibrium at 
around 2 x 10 13 GeV, ?7 S is already suppressed by e~ T ^ T ~ e~ 15 , which justifies the 
approximation of taking the relaxed value rj s = as we did in the two previous sections. 



6. Supersymmetric models 

Baryogenesis and leptogenesis have been also studied in the context of supersymmetric 
theories. The supersymmetry breaking terms decouple at a certain temperature 
T ss ~ 10 7 GeV and supersymmetry becomes a good symmetry above T ss . As a result, 
other global symmetries become effective at high temperatures namely, in the case of the 
minimal supersymmetric standard model (MSSM), the Peccei-Quinn, U(1)a, and U(1)r 
symmetries. It was shown [24J that in addition to B — L another linear combination, R', 
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-3 


1 


-7 


R 


1 








1 

2 


1 

2 


1 

2 


1 

2 


l 

2 


1 

2 





-2 


-10 


Rm 


3 








1 

3 


8 
3 


8 
3 


-3 














-36 



of A, R, B and L is anomaly free under the gauge groups SU(3) and SU(2). Thus, R' 
and B — L are unconstrained conserved quantum numbers even if all reactions including 
sphaleron processes are in thermal equilibrium. Then, a net baryon number asymmetry 
requires non-zero values of either R' or B — L but not necessarily both of them, in 
contrast with the non-supersymmetric case where a non-zero B — L is needed. 

But at very high temperatures there are in fact other approximate symmetries 
because some of the lepton and quark Yukawa couplings are thermally decoupled 
as discussed in the previous sections. In addition, as we prove in the following, 
supersymmetry adds two and not just one SU(3) and SU(2) anomaly free conserved 
quantum numbers. 

The exact supersymmetric part of MSSM comprises the SU(3)xSU(2)xU(l)y 
gauge interactions and the Yukawa couplings specified by the superpotential 

W = \ u H 2 qu c + \ d H iq d c + A e #i/e c , (36) 

where the flavour indices are omitted. In the case of supersymmetric leptogenesis 
123 1213 I2H I2H] one adds heavy sterile neutrinos N a with Majorana masses M a and the 
superpotential 

W N = ^MNN + X N H 2 1 N . (37) 

The MSSM interactions conserve B and L and two other quantum numbers, A and R, 
of the global Peccei-Quinn and U(l)# symmetries. The particle quantum numbers are 
specified in table ^(i? as defined in (2E1I2II])- Hi and H 2 are the Higgsino fields and A 
stand for the left-handed components of the gauginos W, B and g. 

The R charges of the Higgs fields H\, H 2 and scalar superpartners of q, u c , d c , 
I, e , N exceed the respective fermion superpartner charges by 1 unit so that the 
superpotential carries 2 units. The gauge bosons W, B and g have zero A and R 
charges. Rm is the linear combination 

R m = 3R+1b-^L, (38) 

defined in such a way as to be non-anomalous under SU(3) and SU(2), and conserved by 
sterile neutrino Majorana masses and Yukawa couplings, if applicable. The Rm charges 
of the Higgs and matter scalar particles exceed their fermion superpartner charges by 3 
units and the superpotential carries 6 Rm units. 
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The last three columns of table Q represent the anomalies under the standard 
model gauge groups SU(3), SU(2) and U(l)y, respectively, normalized to coincide 
with the quantum number variations per unit of Chern-Simons number, AN. For 
example, the baryon number varies as AB = 3AN under SU(2) weak sphalerons while 
AA = AN. The Higgsinos and gauginos contribute to the anomalies as follows: the 
Higgsino fermion numbers vary as ANg = AN g 2 = AN under weak sphalerons and 
the Wino left-handed components as AN^y = 4AN. The gluino left-handed components 
vary as ANg = 6 AN under SU(3) QCD sphalerons. As a result the chemical potential 
constraints enforced by weak and strong sphaleron processes when they are in thermal 
equilibrium are modified with respect to the standard model. Equations (jHJ) and (pTTj) 
are replaced with 

Vw = ( 3r ^ + Vh) + m t + Vh 2 + ^Vw > ( 39 ) 

i 

Vs = Yl ( 2? fe ~ ^ ~Vdi) + 6r lg ■ (40) 

% 

In order to determine the rj w and rj s damping rates it is necessary to identify 
the quantum numbers conserved by the respective sphaleron processes. Besides the 
hypercharge and B — L there are flavour dependent quantum numbers that are conserved 
by both weak and strong sphalerons. Within the standard model the partial lepton 
numbers of the lepton iso-singlets, L Si , the differences between the right-handed charm 
or down quarks baryon numbers and the Ur baryon number, B CR —B UR and B d .—B UR , the 
differences between total lepton and baryon numbers of different generations, Li — Lj and 
Bi — Bj are only violated by Yukawa couplings. At high enough temperatures these are 
out of equilibrium and the respective flavour quantum numbers are conserved to a very 
good approximation. This remains true in the supersymmetric version of the standard 
model simply by extending the partial lepton and baryon numbers to the s-leptons and 
s-quarks scalar superpartners. But there are two more conserved charges. One is Rm 
as discussed above. The Peccei-Quinn charge A, on the contrary, is not conserved by 
strong sphalerons, nor any combination of A, B and L. However, a combination like 

A u = A- l -B- 9B UR (41) 

is not anomalous under SU(3) and SU(2) and is conserved as long as the up quark 
Yukawa coupling remains out of equilibrium, which happens at temperatures above 
10 8 GeV. Thus, supersymmetry adds two new conserved charges, A u and Rm, at very 
high temperatures. 

Another significant consequence of supersymmetry is a different evolution of the 
gauge coupling constants which reflects on the sphaleron rates that scale with the fifth 
powers of ct s and a w . At 10 12 GeV, a w ~ 1/26 instead of the standard model value 
a w ~ 1/41. It makes the weak sphaleron diffusion rate 10 times larger than in the 
non-supersymmetric case: r w rs 25 a% T 4 w 2 x 10~ 6 T 4 . A similar increase occurs for 
QCD sphalerons. At 10 14 GeV, a s « 1/21.5 yields a diffusion rate, T s ~ 300a s 5 T 4 w 
6 x 10~ 5 T 4 , 15 times larger than in standard model case (a s ~ 1/39). 
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6.1. QCD sphalerons 

The determination of the strong sphaleron degeneracy parameter r] s proceeds like in 
the standard model case. One writes down the chemical potentials in terms of a set 
of independent global quantum numbers that are conserved by all significant processes 
with the possible exception of QCD sphalerons and obtains a relation like (jHS| for the 
quantity defined in equation P0|): 

Vs T 3 V = -aB UR + b a Q a . (42) 

Q a are all charges conserved in a temperature range around 10 14 GeV. If the MSSM 
interactions are the only ones in thermal equilibrium, the charges Q a are the baryon 
number B, only marginally violated by weak sphalerons, and the anomaly free B — L, 
Rm, A u , L e ., B CR — B UR , B di — B UR , Li — Lj and Bi — Bj. If other non-standard 
interactions are in equilibrium new constraints apply and the number of conserved 
charges is reduced. In any case the variation of i] s T 3 V is always the one of B UR times 
the factor —a. Equation (|34j) and (J3*5j) still hold and the relaxation temperature is given 
by T* = aT s /6. What are different are the coefficients a and b a , which depend on the 
actual degrees of freedom and interactions in equilibrium. 

The gauge and gaugino interactions put constraints of the general form 

= Ml + Mi , (43) 

where fi stands for any of the left-handed fermion particles q i} uf , df , l iy ef, Hi, H 2 
and fi is the respective scalar superpartner. = = ^ — ^ is the common 
(left-handed) gaugino chemical potential whereas the gauge bosons have zero chemical 
potentials. The other constraints are enforced by the top quark Yukawa coupling, 

/V - Vt R +Hh 2 = 0, (44) 

and the condition of zero total hypercharge. In the usual linear approximation this 
amounts to: 

+ 2 /V - Vdi - - Me J + Vh 2 ~ Phi 

i 

+ 2 + 2 M«, - - ~ P*) + 2 ^h 2 -n Hl ) = 0- (45) 

i 

The constraint equations can be used to express all chemical potentials in terms of 
B UR and the conserved quantities Q a of equation (JUJ). If no interactions beyond 
the MSSM are in equilibrium one obtains the coefficient a = 9440/429 and the 
relaxation temperature T* = 3.67 T s . The temperature T s = r s T~ 2 i/ -1 calculated 
with T s ~ 300 a s 5 T 4 , a s ~ 1/21.5, and an effective number of degrees of freedom 
g* ta 229 (H = 1.66 gl /2 T 2 / 'M Pi ) amounts to T s ~ 3 x 10 13 GeV. It makes the relaxation 
temperature, T* ~ 10 14 GeV, 4 times larger than in the non-supersymmetric case. 
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6.2. Weak sphalerons 

When the weak sphalerons approach thermal equilibrium the strong sphaleron processes 
are already very rapid and the respective degeneracy parameter is damped down by 
several orders of magnitude. One applies then the constraint ?7 S = on the expression 
(140 j) and the constraints (f4*3*|) and (|44jl . The particle chemical potentials and the weak 
sphaleron degeneracy parameter r] w of equation (|3U|) can be written in terms of the 
baryon number and a set of quantum numbers Q a conserved by both weak and strong 
sphalerons: 

Vw T 3 V = aB + b a Q a . (46) 

If no other than MSSM interactions are in thermal equilibrium the Q a charges are 
B — L, R M , A u , L ei , B CR — B UR , B d . — B UR , L { — Lj and Bi — By As a first estimate 
of the ?7 W relaxation temperature one ignores the violation of L TR and Bb R — B UR by the 
tau and quark Yukawa couplings. Then, the time derivative of i] w T 3 V is just adB/dt 
and the results (fTTj) and (|T2"j) follow from there. We obtain a = 8179/4720 and the 
relaxation temperature T* = |aT w = 2.60 T w . Taking a diffusion rate T w « 26a^T 4 
with a w » 1/26, one gets T w = V^T^H' 1 « 10 12 GeV and T* » 3 x 10 12 GeV. 
Although T w is 7 times larger than in the standard model, the relaxation temperature 
T* is just less than 2 times larger. It remains so in the vicinity of the 6r and tr Yukawa 
coupling relaxation temperatures causing the same phenomenon discussed in section 0] 
of damped oscillations of rj w and the other two degeneracy parameters. 

7. Conclusions 

Weak and strong sphaleron processes enforce chemical potential constraints at high 
temperatures but the respective relaxation temperatures are not the ones where the 
sphaleron frequencies TT~ 3 coincide with the Hubble expansion rate. They are larger 
than that for both kind of sphalerons but the exact factors depend on the set of reactions 
in thermal equilibrium at the relevant epoch: they are 10 times larger in the case of the 
standard model and around 3 in the case of supersymmetric models. These numbers do 
not depend on the origin of the asymmetries i.e., on the particular lepto-baryogenesis 
mechanism, except if exotic particles or interactions are in thermal equilibrium at the 
time considered. 

One aspect that has been overlooked so far concerns the evaluation of the sphaleron 
diffusion rates themselves. The a w and a s gauge coupling constants run with the energy 
scale and this causes the weak and strong sphaleron rates to approach each other at very 
high temperatures and to be model dependent on the other hand. The weak sphaleron 
rate is about 10 times larger in supersymmetric model extensions than in the standard 
model while the strong sphaleron rate is enhanced by a factor of 15. Our calculations 
yield the relaxation temperature values of 2 x 10 12 GeV and 3 x 10 13 GeV for weak 
and strong sphalerons, respectively, in the standard model case and 3 x 10 12 GeV and 
10 14 GeV in the supersymmetric version. 
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Weak sphalerons enter in equilibrium in the same period as the bottom and tau 
Yukawa couplings. We shown that the respective degeneracy parameters obey a coupled 
system of equations and oscillate into each other as they are damped away. 
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